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Abstract: We consider infinite dimensional port-Hamiltonian systems. Based on a power 
balance relation we introduce the port-Hamiltonian system representation where we pay 
attention to two different scenarios, namely the non-differential operator case and the differential 
operator case regarding the structural mapping, the dissipation mapping and the in/output 
mapping. In contrast to the well-known representation on the basis of the underlying Stokcs- 
Dirac structure our approach is not necessarily based on using energy- variables which leads to 
a different port-Hamiltonian representation of the analyzed partial differential equations. 
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1. INTRODUCTION 



Modeling of physical systems described by partial differ- 
ential equations (pdes) in a port-Hamiltonian framework 
has been treated extensively over the last years. However, 
there are still many open issues and there exists no unique 
representation. An important contribution in this field was 
the introduction of Stokes-Dirac structures which allows 
to analyze field theories in a port-Hamiltonian framework 
and to exploit this system representation for the controller 
design, see e.g. van der Schaft and Maschke (2002); Mac- 
chelli et al. (2004a,b); Macchelli and Melchiorri (2004); 
Maschke and van der Schaft (2005). Roughly speaking, the 
key property of the Stokes-Dirac structure is to represent 
the power balance relation of physical systems in an formal 
way. This is achieved by combining so-called flows and 
efforts in the domain and on the boundary of a system. 
Using this approach the proper choice of energy variables is 
crucial and the Hamiltonian is also a functional depending 
on energy variables. 

A different kind of a port-Hamiltonian representation, 
not focusing on energy variables, but also based on a 
power conservation law, was proposed in Schlacher (2008); 
Schoberl et al. (2008). The proposed Hamiltonian frame- 
work can be seen as an extension of Olver (1986) by incor- 
porating dissipation and boundary ports. This formulation 
has been exploited for control issues in Schoberl and Siuka 
(2011) and Siuka et al. (2011). 

In this contribution we will present an enhancement of 
Schlacher (2008); Schoberl et al. (2008) in such a way that 
beside the input map as in Schoberl et al. (2008) also the 
structural mapping J and the dissipation mapping 1Z may 
involve differential operators. 



In the approach based on Stokes-Dirac structures it is 
quite natural that these mappings involve differential 
operators, however due to the choice of energy variables 
the Hamiltonian density does not depend on derivative 
variables, in contrast to the approach presented in this 
paper. For instance in mechanics, the choice of energy 
variables suggests to use the strain e as the dependent 
coordinate, whereas also the displacement u can be used. 
This latter choice leads to derivative coordinates when the 
potential energy has to be stated, since the energy is a 
function of the strain and e = dxu is met in the case of 
a one-dimensional domain. This has severe consequences 
in the application of the variational derivative applied 
to the Hamiltonian density, leading to a different port- 
Hamiltonian representation, and has also an impact on 
the computation of structural invariants, see also Schoberl 
and Siuka (2011). 

This paper is organized as follows. In Section 2 some 
notation is presented and the geometric objects that play 
a fundamental role in the paper are introduced. The 
third section is dealing with the representation of port- 
Hamiltonian systems described by pdes, where we focus 
on two cases: the non-differential operator case and the 
differential operator case. Two specific applications, the 
vibrating string and a simple model of magnetohydrody- 
namics are analyzed in the fourth section to demonstrate 
how the introduced differential operators have an impact 
on the power balance in practice. The paper closes with 
some concluding remarks. 

2. NOTATION AND PRELIMINARIES 

In this paper we will apply differential geometric methods 
and we will use a notation that is similar to the one 



in Giachctta ct al. (1997). To keep the formulas short 
and readable we will use tensor notation and especially 
Einsteins convention on sums. We use the standard symbol 
® for the tensor product, A denotes the exterior product 
(wedge product), d is the exterior derivative, J the natural 
contraction between tensor fields. By arc meant the 
partial derivatives with respect to coordinates with the in- 
dices B . Furthermore C°°(-) denotes the set of the smooth 
functions on the corresponding manifold. Moreover we will 
not indicate the range of the used indices when they are 
clear from the context. Additionally pull back bundles are 
not indicated to avoid exaggerated notation. 

In the forthcoming wc will consider bundle structures in or- 
der to be able to separate dependent and independent co- 
ordinates. Let us consider the bundle X — > V, (X A , x a ) — > 
(X A ) where x are the dependent and A the independent 
coordinates. The first jet manifold J X {X) can be intro- 
duced possessing the coordinates (X A ,x a ,x A ), where the 
capital Latin indices A, B are used for the base manifold 
V (independent coordinates) and x A denote derivative 
coordinates of first order (derivatives of the dependent 
coordinates with respect to the independent ones). The 
jet structure also induces the so-called total derivative 

d A =d A + x A d a + x AB d® 

acting on elements including first order derivatives and 
x AB correspond to derivative coordinates of second order 
living in J 2 (X), the second jet manifold. We will mainly 
focus on the first order case in the sequel, however the jet- 
structure and the total derivatives can easily be extended 
to higher order cases. 

We will treat so-called densities in the sequel (a quantity 
that can be integrated), where we pay special attendance 
to densities of the form £ = TdX with T E C°°(J rl (A')) 
where dX denotes the volume element on the manifold 
V, i.e. dX = dX 1 A ...AdX d with dim(D) = d. Since 
T 6 C°°{J 1 {X)) is met in our case we restrict ourselves to 
first order theory and additionally we denote by F — f v $ 
the integrated quantity, where of course the map x = <&(A) 
leading to xa = 8a^(X) has to be plugged in to be able 
to evaluate the integral properly. 

Based on the bundle structure X — > V one can introduce 
several tangent structures, where for our considerations 
the vertical tangent bundle V(X) and the cotangent bun- 
dles representing differential forms will be used. Of special 
importance will be the space Af (X) induced by the bundle 
structure X — > V 

Af(X)=T*(X) A (ATOP)) 

with a typical element u> = uj a dx a A dA. It is worth 
stressing that the functions uj a may depend on derivative 
coordinates, however as stated above to simplify the nota- 
tion no pull backs will be indicated in the forthcoming, 
i.e. it should be clear from the context which order of 
derivative is included in the differential form 1 . An im- 
portant object is the horizontal exterior derivative d/j, 
which meets dh(4>) = dX A A d A {4>) acting on a differential 
form 4>, where d A {(f>) denotes the Lie-derivative of 4> with 

1 To be more precise: If u) a e J X (X) then u> e (n^)* Af(X) with 
7Tq : J' 1 (X) — > X. For simplicity we write u) g Af(X) and suppress 
the pull back. 



respect to d A (see the appendix for more details concerning 
the relationship of d and dh and Stokes theorem in that 
context). 

Furthermore we will treat linear differential operators (of 
order k) that are of the following form 

T> : Af(X) -> V(X) 

that maps an element Af (X) of jet-order p to an element 
V(X) of jet-order p + k. In coordinates we have 

3>(w) = ^ a/3M d A (w a )d , d A = d Ak o . . . o d Al . 

with lj G Af(X). The adjoint operator £>* follows by 
integration by parts and fulfills the condition 

D{lo)\vj = T)*{zu)\lu + d h {d) (1) 
with iv, w € Af(X), where is a bilinear expression 
involving the total derivatives up to order k — 1, see Olver 
(1986). 

3. SYSTEM REPRESENTATION 

In this section we will introduce port-Hamiltonian systems 
described by pdes based on a power balance relation. This 
means that the system is introduced in such a way that 
the power balance relation together with the structure 
of the equations represent the physical process. We will 
consider two cases: (i) the Hamiltonian depends on deriva- 
tive coordinates, but the operators J 1Z and Q (namely 
the interconnection, the dissipation and the input/output 
maps) are just linear maps and (ii) where we relax the 
assumption regarding J, 1Z and Q and consider adequate 
differential operators 3, 9^ and 25 instead. Before we are 
able to introduce the corresponding system representations 
we need to introduce some geometric concepts which will 
be exploited later on with respect of the derivation of the 
power balance relation. 

Let us consider a vector field which is used to measure the 
change of the density T . We use a (generalized) vertical 
vector field v : X — > V(X) locally given as v — v a d a where 
v a may depend on derivative coordinates, together with 
its first jet- prolongation ^(v) which reads as 

j 1 (v)=v a d a +d A (v a )d£. 

Then we compute the Lie-derivative of the densitiy J as 
it has been defined before with respect to the vector field 
^(v) and we obtain the important relation 

j\v){TdX) - (v a (d a F - d A d*F) + d A (v a d£F)) dX 

= (v a 5 a F + d A {v a d A F)) dA. (2) 

Here the Euler Lagrange operator 5 

5 : J 2 (X) -»• Af (X) (3) 

comes into play, whose coordinate expressions involves the 
variational derivative S a which acts on T as 

5 a F = d a T - d A d A T. 

Remark 1. In fact (3) written out in coordinates is a map 

5 a F — > 5 a Fdx a A dA 

with 3 = TdX. 

Applying the Theorem of Stokes (see again the appendix) 
to (2) we find that 



f j 1 (v){TdX) = [ v a (S a T) dX+ [ v a (d^T) dX A 

Jv Jv JdV 

= [ v\S$+ [ v\5 9 $ (4) 

Jv JdV 

with dXA = cttJdAT and the boundary operator 

5 a $ = d£j r dx a AdX A . 
The relation (4) will be of key interest in the forthcoming, 
since it provides a natural decomposition of the expression 
Jjy j 1 (i>)(J : cLY) into a term on the domain V and one the 
boundary dT>. 

Important is furthermore the case when the generalized 
vector-field v is linked to the solution of a pde system 
(via its semi-group, that v may generate), then the formal 
change of F = f v TdX along solutions of a pde system can 
be computed as J T) j 1 (v)(FdX) (provided all operations 
are admissible), which will be denoted by 

F = [ j\v){TdX) 
Jv 

in this special case. 

3.1 The non- differential Operator case: 



We introduce a port controlled Hamiltonian system on a 
bundle X — > T> in the non-differential operator case in the 
form of 

x = (J — 1Z)(5$j) +u\Q 

v = g*m (b) 

together with appropriate boundary conditions and with a 
first order Hamiltonian Sj = HdX with H G C°°( l 7 1 (A')). 
(Here it is worth mentioning that we restrict ourselves to 
the first order case). 

Remark 2. It is worth stressing again, that of course 5 
is a differential operator, and the phrase 'non-differential 
operator case' only refers to the maps J, 1Z and Q. 

We make use of the linear maps (over vector bundles) 

J,K: Ai^V(X) 
where J is skew-symmetric and meets 

J(oj)\(m)+J(vj)\(Lu)=0 (6) 

and since 7Z is symmetric and positive semi-definite we 
have 

TZ{lu)\(w) -U(zu)\(lu) =0, K(lo)\(lo)>0 (7) 
with uj,-cu e Af . 

The input and the output bundles are given as U — > X 
and y — > X and the map Q and its adjoint Q* are of the 
form 

g-.u^v(x), g*-.Af^y. 

with 

{ u \g)\oj = u\(g*\uj). (8) 

Replacing $ by f) in (4) and setting v — x we obtain 
H= [ (J-Tl)(SSj)\SSj+ [ (u\g)\6f)+ [ x\5% 

JV Jv JdV 

= - f K(6Si)\6Si+ f u\y+ [ x\S% (9) 

Jv Jv JdV 



which reflects the power balance, since the total change 
of the functional H along solutions of (5), is made up of 
dissipation, collocation on the domain and an expression 
corresponding to a boundary port (if it exists) depending 
on the boundary conditions, see Schobcrl ct al. (2008); 
Siuka (2011). 

Remark 3. Since j v u\ y in coordinates reads as JjjtfyidX 
we interpret the outputs as density valued such that the 
fibers of the vector bundle U — > X possess the base e, and 
the fibers of the vector bundle y — >• X possess by duality 
the base e l ®dX. Similarly, this is reflected also in the case 
of the pairing of V{X) with Af which also takes its values 
in the densities since with v £ V(X) and co S Af we have 
v\uj = v a uj a dX. 

Casimir junctionals Let us consider how to obtain so- 
called Casimir densities (or functionals) where we restrict 
ourselves to the first order case only, i.e. € = CdX with 
C G C°°(J rl (A')). Using (4) where we replace T by C, 
setting v = x and C = f v £ the relation 



C= / (J-K)(5f))\5£+ / (u\g)\6£- 

Jv Jv JdV 

is obtained and a Casimir functional meets 



x\S 9 € 



C= / (u\g)\5€. 
Jv 

This leads to the following two conditions for the Casimir 
density 



5 a C{J af) -U af} ) = Q 



(x a d£C)\ 



dv ' 







(10) 

(11) 



which have to be fulfilled (independently of the Hamilto- 
nian density Sj). If in addition 

/ (u\g)\S€ = o 

Jv 

is met, then the density is a conserved quantity. 

The concept of control by interconnection based on 
Casimir functionals using the presented system representa- 
tion has been applied to a heavy chain system in Schdberl 
and Siuka (2011) and to the Timoshenko beam in Siuka 
et al. (2011). The generation of Casimir functionals based 
on the Stokes-Dirac structure and its consequences for 
control have been treated for instance in Macchelli et al. 
(2004b); Pasumarthy and van der Schaft (2007). 

3.2 The differential operator case: 

A port-Hamiltonian system with respect to the differential 
operator case reads as 

x = (3 - 3t) {SSj) + 25 (u) (12) 
y = (&*(SS)) 

together with appropriate boundary conditions, see also 
Siuka (2011). 

For this case the operators 3,91 :Af — >■ V(X) are r-order 
linear vector valued differential operators where 3 is a 
skew-adjoint operator according to 

3(w)Ja7 + 3(a7)Jw = d fc (j) (13) 

with 



)=) A dX A , uj,weAf, 

cf. (1). Furthermore, 91 is a non-negative self-adjoint 
operator, i.e., 

m(Lu)\zu - 9t(ro)Jw = d h (r) , 9t(w)Ju;>0. (14) 
with c = v A dXA- 

The input operator © and its adjoint operator 25* are maps 
according to 

<5:U^V(X), &*:Af^y (15) 
and they are linear r-order differential operators with 
respect to the relation 

0(u)Jw = uJC*(w) + d fc (fl) , Q = Q A dX A , (16) 
Consequently, from (16) we are able to derive the relation 

<5{u)]5Sj = u\e* (5Sj) + d h ( Q ) 

= u\y + d h ( Q ) (17) 
where the pairing u\y will correspond to the domain port. 
Remark 4- It should be noted that the relations (13), (14) 
and (16) are the generalizations of (6), (7) and (8) to the 
differential operator case, where especially the horizontal 
derivative will lead to appropriate boundary terms when 
the adjoint operator is taken into account. 

The formal change of the Hamiltonian functional H along 
(12) now takes the form 

H= [ (3-9t)(J$)J<H5+ / e(u)\SSj+ [ x\S d Sj (18) 

JT> JV JdV 

and due to the involved differential operators it cannot be 
concluded in a unique fashion how the power is extracted 
over the domain and/or the boundary. It is easily checked 
for instance that plugging in (17) into (18) leads to colloca- 
tion on the domain and an additional boundary expression. 
These non-uniqueness depends on the properties of the 
used differential operators in 3, 9^ and 25 which also has 
an impact if Casimir functionals are to be considered in 
this differential-operator case. To be able to draw some 
conclusions also in that case we will restrict ourselves to 
specific operators described in the following section: 

Specific Operators Motivated by the forthcoming appli- 
cations we will consider two types of operators. 

We introduce a second-order non-negative self-adjoint op- 
erator 91 locally given as 

m(iu) = d A (m^ AB d B (ujp))d a , (19) 

with 

g^afiAB fj^fSaBA <y^af3AB g (J°° {X) 

It satisfies the relation 

x A dX A = m(uj)\zu d -m{zu)\uj d , (20) 

cf. (14) with 

Wd = -d A \u, wg = -d A \w 
where we used 

9t(w)=9t a ' MB d B (oj p )d a . 

Remark 5. The condition <R a P AB G C°° (X) might be 
relaxed, and C°° (X) can be replaced by an appropriate 
jet-space. Since, in our concrete application this is not the 
case, we will not focus on that fact more detailed. 



The operator 9^ is non- negative since 91(cj)Jcj can be 
written as 

-d A {u a )^ AB d B {up)dX + d h {^) (21) 

with the boundary expression 

<& m =u a V\ afiAB d B (u p )dX A . 
The non-negativity of the operator follows if 

-d A (u a )y\ afsAB d B {up) >0 (22) 

is met. 

Furthermore, we define a first-order input operator © with 
25 (u) = <8? A d A (u l ) d a (23) 
and, <&f A e C^iJ^X)) which meets 

<8{u)\cj = u\e*(i J j) + d h (e(u)\LU d ) (24) 

with 

(u) = <3* A u l d a , uj d = -d A \uj. 
The relation (24) takes the form of 



co a <3f A d A (u t )dX = -u*d A {u; a e? A )dX + d, l (» )(25) 



uA\ 



with 



=w a <8f A u l dX A . 

Remark 6. As above the condition & f A e J X {X) can be 
relaxed in the same spirit as for the dissipation operator. 

4. EXAMPLES 

In this section we will discuss the presented ideas based 
on physical examples. We will present the vibrating string 
example in two fashions -with and without structural 
damping which affects the 9i operator. To show how an 
input operator influences the energy flows we analyze a 
simplified version of magnetohydrodynamics, similar as we 
did in Schoberl et al. (2010). 

4-1 Vibrating string 

Let us consider the system of a vibrating string. As 
independent coordinate we choose the spatial coordinate 
X, the dependent coordinates are the deflection w and 
the temporal momentum p. This leads to the following 
bundle X — >• T>, (X,w,p) — > X, The first jet manifold 
J 1 (X) additionally includes the derivative coordinates wx 
and px and the boundary dV consists of two points only, 
namely X = and X = L where L is the length of the 
string. Approximately, the pdes take the form 

" = P (26) 
p = dx(P(X)w x ) 

with additional boundary conditions and the force in the 
string is denoted by P(X). To rewrite this system in a 
Hamiltonian fashion we consider the Hamiltonian density 
fi = HdX with 



H 



(27) 



The total energy can be evaluated by H = J HdX. To 
obtain partial differential equations in the form as in (5) 



we can set 1Z and Q to zero (no damping and no inputs 
acting on the domain). Then we easily obtain 



(28) 



w 




' o r 




~S W H~ 


p _ 




-1 




5 P U_ 



J(SSj) 

with x = (w,p) together with the expression 

s w n = (d w - d x d*) n = -d x (p(x)w x ) 

and 

S p H=(d p -d x d p x )H = ^. 
Evaluating (9) gives 



H= [ x 

JdV 



x]S 9 Sj= wP(X)w x \^ 

and again the power balance relation where the boundary 
conditions are crucial in the determination of that expres- 
sion. 

Now we extend the equations by structural damping r > 

as 



. P 
w = — 



P^ 



w 




' o r 




p _ 


=( 


-1 







~ 5 W H~ 


]) 


6 P H 



p = dx(P(X)w x ) + d x (rd x -). 

P 

The port-Hamiltonian representation can be given as 

'0 

-d x (rd x (-)) 

where now a differential operator is utilized in the d\ 
mapping, as introduced in (19). The power balance reads 
in this case as 

H = - [ V\(5$j)\Sf)+ [ x\5 d S). 

Jv JdV 

Using (21) we obtain 



/ m(6f))\5$j = [ wd x (rd x (w))dX 
Jv Jv 

where the right hand side of (29) decomposes into 



(29) 



/ d x (w) r d x (w) dX ■ 
Jv 



av 



■ d x (w) dX A 



and thus 

H = - d x (w) rd x (w)dX + I <B S 

Jv JdV 

with the boundary expression 

<8 S = w (r d x (ib) + P(X) w x ) dX A . 

It is obvious that the structural damping extracts power 
as can be seen from the term on the domain 

— I dx (w) r dx (w) dX = -r / (d x (w)fdX 
Jv Jv 

in H since r > 0. 

4-2 Magnetohydrodynamics (MHD) 

In this section we consider an example that enables us to 
introduce a differential operator in the in/output mapping, 
namely so-called inductionless MHD. Roughly speaking, 
we analyze the macroscopic behavior of an electrically 



conducting fluid in the presence of external electromag- 
netic fields. We assume that the dynamic of the addition- 
ally induced electromagnetic parts can be neglected (low 
magnetic Reynold's number). Furthermore, we confine 
ourselves to the case of negligible viscosity and electrical 
conductivity. 

To obtain the mathematical model we consider again a 
bundle structure of the form #->•£>, (X A ,q a ) -> (X A ) 
and we introduce the so-called material mass and charge 
density, p{X) and n{X). Furthermore we assume the 
existence of a stored energy function E st (g^) which meets 



P = -p 



dE st 
dJ 



J = det (F) 



with the deformation gradient F% = q% and F^F^ = 6%, 
where we have introduced the material pressure V (X). 
The electrostatic potential is denoted by Aq (q) and the 
vector potential by A a (q), where the components of the 
electric field strength are given by E 0a = d a A and of the 
magnetic flux density by 

B a p = d a Ap - dpA a . 
The Hamiltonian density f) = HdX follows as 

n = ^-6 ap (p a - pA a )( Pfj - fiAp) + P E st 
zp 

with the momenta 

V<x= f>5 a pq p + p.A a . 
Remark 7. The Hamiltonian can be derived using the 
Lagrangian 

£ - \p#*& a piP - PEst +P(A + q a A a ) 

together with H = (p a q a — C) where the potential Aq is 
neglected in a first step, since it serves as the input, which 
is considered separately. 

The partial differential equations read as 

p a - -JF A d A V + p- p 5 pp ( Pp - pA p )d a A p 

and its Hamiltonian representation follows as 



~q a ~ 




' 5f 




's p w 


Pa 




-€ o_ 







with 5 fi n = d^H = £-U. 

Taking the system input Ao into account we extend the 
system and obtain the following representation 



J 



6?H 



and 



© (A ) = 



+ © (A ) 







pF^d B A 



where the input map corresponds to the definition of 
equation (23) with u = A . 

The power balance then again takes the form 
H= [ <5(Ao)\5Si+ [ x\5% 

Jv JdV 

where the interesting term corresponds to the pairing on 
the domain including the input and we obtain 



f <5(A )\6Sj= f ^(d B A )q a dX 
Jv Jv 

= - [ d B (F*S a )A dX+ f Q5 M 

Jv JdV 

with 

< B M = S a F^A dX B , 

see also (25). Here the adjoint operator on the domain 
takes the form 0* (5Sj) = -d B (F^S a ). 

The expression 

S a = ^ {pp - nAp) 
P 

is equivalent to the components of the convective current 
density which is caused by the convective transport of 
charge. 

5. CONCLUSION 

Based on the relation (4) that describes the power balance 
relation, we introduced port-Hamiltonian systems modeled 
by partial differential equations, in such a way that energy 
conservation, dissipation, collocation and energy ports 
come apparent. This key idea is very similar of that used 
in the definition of the Stokes-Dirac structure, however 
due to a different choice of coordinates (not necessarily 
energy coordinates) the port-Hamiltonian representations 
of the two approaches differ significantly. This can be 
observed by a close inspection of the differential operators 
that are needed to describe the physical systems properly. 
Therefore it is of great interest for future investigations 
two analyze and compare these two promising directions 
more closely. 

APPENDIX 

Les us consider a bundle n : X -> T>, (X A ,x a ) -> (X A ) 
and its n-th order jet manifold J n ( X) as well as the bundle 
7Tq : J n (X) -> X together with a section 4> : V -> X, i.e. 

x = 4>(X). 

The exterior derivative d is connected to the horizontal 
derivative dh through 

(j n+ V)*(d h ( W ))=d((j n 0)*M) (30) 

for a form ui living on J n (X). Roughly speeking the pull 
back (j n (f>)* (co) denotes cj o {j n <p) where j n cj> corresponds 
to the n-th order jet-lift of the section <j>. 

Furthermore we have 

/ r +i {4>r{A h u)= f d(r(4>y»)= / r (<«» 

Jv Jv JdV 

for lu e (n r ^)*(A d ^ 1 (X)), which is nothing else than Stokes' 
theorem adapted to bundles. 
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